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1. Introduction 

In the construction of complex modular simulation 
systems, as the one used in High Energy Physics, the 
question of dividing the system to functional parts is 
essential. In the modular approach, the problem can 
be divided into parts, and each part can be addressed 
by different researchers or teams. However, such ap- 
proach is only possible if the mathematical structure 
of the problem have certain algebraic properties. 

In practice, the module-based work model is often 
an idealization: creating scientific software can not be 
separated from the research itself. As a consequence, 
the architecture of programs needs to be modified at 
various steps of projects development, for example to 
accommodate more precise models. 

Another class of difficulties in the development of 
a simulation segment arise from the constraints, im- 
posed by other segments of the large project, namely 
in the definition of data structures and interfaces. 

Finally, the demands of the end-users, their exper- 
tise in handling simulation blocks (which is often lim- 
ited) and understanding the conceptual models used 
are of critical importance. 

These types of difficulties are quite universal for 



any complex, scientific problem. In the following let 
us concentrate on a relatively simple model (yet al- 
ready quite sophisticated) for the simulation of QED 
radiative corrections in decays, and the corresponding 
simulation package PHOTOS II 11213 141 . We believe that 
presented "development drama" can be of interest not 
only to the readers interested in QED bremsstrahlung, 
but also in general case. In this respect our presenta- 
tion can be understood as a summary and invitation to 
reading Q. With respect to previous |6| longer ver- 
sion of the presentation let us point on the particular 
property of the simulation of decays and in particular 
necessary in that step radiative corrections. Function- 
ality of such modules in some cases can be attributed 
to the physics of hard processes but in the other ones 
to the periferic ones, often to be used to measure spin 
states of the decaying particles. In the second case 
lower energy experiments provide invaluable frame- 
work necessary to evaluate systematic error for high 
energy applications. Finally such simulation modules 
need to match (directly or with the help of tests) the 
ones of low energy e + e~ experiments. 

Our contribution is organized as follows. We start 
with the presentation of the PHOTOS algorithm in sec- 
tion 2. We highlight the role of the event record in 
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program's construction. We will also specify those 
properties of event record that are necessary for the 
program precision. 

Section 3 is devoted to the description of recent 
years' developments in the PHOTOS algorithm, which 
lead to the improvement of its performance, but at 
the same time introduced more strict constraint on the 
event record used as a data source. We will present 
our efforts in preventing the complications for other 
users of data structures. 

During the years of project's evolution special tech- 
niques devoted to detect and resolve various types of 
difficulties related to event record(s) were developed. 
Section 4 will be devoted to their presentation. In- 
ternal instrumentation utilities of PHOTOS: debugging 
subroutine PHLUPA and kinematic rounding error cor- 
rection subroutine PHCORK will be discusses. Then, 
MC-TESTER [7| will be briefly presented. Finally, we 
will justify the need to duplicate the event information 
which is originally stored in a standard event record 
HEPEVT. 

Section 5 will be devoted to the discussion of the 
issues related to event record standardization. The 
question of porting the algorithm to C++ - the process 
technically completed in 1999 - [8| will be covered 
in that section as well. We will stress the important, 
yet often ignored, aspect of event record construction: 
the structure must not only be convenient and flexible 
from the software engineer point of view, its contents 
must also be clear, from the point of view of involved 
physics models, to the end-users processing the data 
produced by Monte Carlo event generators. This is 
particularly important if extensions of the standards 
are to be not only proposed but also used. 

The summary, Section 6, closes the paper. 

2. Basic design of PHOTOS 

Already at an early step of preparation for the x- 
lepton polarization measurement at LEP1 it became 
evident [9|, that bremsstrahlung corrections in x de- 
cays are necessary for proper modeling of theoretical 
predictions for measurement of the Z boson couplings 
using x-polarization. For that purpose, special rou- 
tine RADCOR was designed. With certain probability it 
was replacing the decay products of x simulated using 
TAUOLA iflOl by the ones with an extra photon added. 
This action was performed in a strict environment of 



the explicit list of momentum four-vectors and in the 
rest-frame of x. In fact, the parameter defining the ac- 
tual decay mode of x was also passed into RADCOR, 
identifying the physical process. 

The x decay products, modified in such way, were 
passed into the TAUOLA interface to KORALZ ITT1 . 
which was the Monte Carlo program for x pair pro- 
duction at LEP1 energies. 

This awkward, yet useful, design was only tailored 
for a single application, and was missing the docu- 
mentation. The limits of its reliability were not ex- 
plored at all. Indeed, the fact that it could properly 
simulate the leading-log parts of the first-order cor- 
rections was enough. The issues of dependence on 
the choice of the gauge were not risen, nor were the 
question of the phase-space coverage. 

However, this exercise provided an important ob- 
servation: to simulate the dominant part of QED cor- 
rections it is sufficient to search through the full event 
structure, identify the branchings corresponding to el- 
ementary decays of some particles, extract the infor- 
mation describing all decay products, and apply a rou- 
tine such as RADCOR. In the early 90's and late 80's 
the HEPEVT event record was providing the complete 
environment that could allow for extension in use of 
RADCOR routine. The first version of PHOTOS |Q] was 
created. Its design relied heavily on the assumption 
that the HEPEVT event record hosts a tree-like struc- 
ture and that all pointers to decay products (daugh- 
ter pointers) and origins (mother pointers) are defined 
and consistent. At each decay splitting the energy- 
momentum conservation was supposed to be fulfilled 
exactly, even if only with a single-precision computer 
arithmetics level. 

The PHOTOS generator gained popularity, and in 
the documentation of its 1994 version @ a multi- 
tude of tests were presented for physical processes for 
which the theoretical predictions were available. The 
corrections for double bremsstrahlung were added as 
well. Further development was stalled due to limited 
interest in improvements, expressed by experimental 
users. 

3. Toward high precision in PHOTOS 

Already in [2] it was found, that if the algorithm 
for the single-photon generation is properly iterated, 
the leading corrections of the double-photon emis- 
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sions can be incorporated as well. In J3] this itera- 
tive solution was extended to multiple-photon emis- 
sion and the tests for leptonic Z decays have shown 
that this solution reproduce the results for the final- 
state bremsstrahlung of the KKMC Monte Carlo pro- 
gram [12 1 with amazing precision. Few years earlier, 
other tests and extensions important for Higgs or W 
decays were introduced also [13]. The good perfor- 
mance of the program was related to other improve- 
ments, introduced at that time and valid for all de- 
cays. For example the implementation became ex- 
act in the soft-photon limit. These improvements, 
however, are of limited importance from the point 
of view of software organization of PHOTOS interface 
with other packages. 

Let us now concentrate on another class of correc- 
tions, discussed in refs. [4 14 1 and corresponding to 
the implementation of the exact, first-order matrix- 
element kernel in bremsstrahlung correction genera- 
tion for Z — * l + r and B — ► K(n)K(n) decays. In the 
latter case, all possible combinations of charges and 
the replacement of n with K were used. 

Let us start with the presentation of the numerical 
results of the test for Z — > decays taken from 

ref. 0. 

For that purpose we have to analyze the plots giving 
largest discrepancies between PHOTOS (run with stan- 
dard options) and KKMC (run with second order matrix 
element and exponentiation). The plots presenting the 
invariant mass of the pair and the invariant mass 
of the two hardest photons for the events with at least 
two photons of energies above 1 GeV in the rest frame 
of the Z are convenient for that purpose. The rates for 
the event samples predicted by the two programs are 
given in the figure's caption. The agreement between 
KKMC and PHOTOS is better than 0.1 % (if calculated 
with respect to the total Z decay rate), yet the differ- 
ences are still visible from the results of simulations 
with 10 8 events. If, as in figure 2, the complete NLO 
kernel is activated in PHOTOS, the differences get re- 
duced by about a factor of 50! This is rather interest- 
ing result of ref. m. 

At this point it is important to realize what is the 
price to pay for such improvement. Certainly, it is not 
the computer time - it remains small and the samples 
of order of 10 9 could easily be simulated overnight 
(PHOTOS is in fact significantly faster than KKMC). To 
answer this question one has to recall the formula for 



the final Monte Carlo weight in PHOTOS. 

Let us write (separated from the phase-space Jaco- 
bians) the explicit form of the real-photon matrix el- 
ement, as used in the standard version of PHOTOS (as 
published in HEJ) for the e+e~ -> Z°/f -> n+p (j) 
process: 

yPHOTOS _ 6' 2 «(1-A) 2 1 fii/i „, \2\ J 
X f - Arts S V^iT^ + ^-Xk) ) \ 



dc B l J(l-cos0+) .;(l+cos0 + ) 
da \ s ' 2 ' 2 
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where :0+ = Z(p+,q+), 0_ = Z(p-,q-), 
©Y = Z(y,/j~) is defined in (/j + ,// - )-pair rest frame. 

For its calculation (with respect to Born cross- 
section) it is enough to know the four momenta of the 
Z and its decay products. In the presented formulae 
we follow the notations from refs. 041151 . This ex- 
pression is to be compared with the exact one, taken 
from ref. ifTBI : 



y C°«(l-A)^ J , 
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The resulting weight is rather simple, and reads: 



WT 2 = 



A<>B 



( s .,m') + ^(s.,\ u ) 



(Hkoag?) /. 3 a \ 
2 



For its calculation the numerical value of the elec- 
troweak couplings of Z to fermions, as well as infor- 
mation on the state from which the Z was produced 
is nonetheless necessary. This seemingly trivial re- 
quirement puts new requirements on the event record: 
the details of the process of the Z productions need 
to be coded in the event record, then correctly deci- 
phered by PHOTOS to calculate the process-dependent 
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weight. From our experience this requirement of 
PHOTOS may be difficult to accept by other users of 
event records. The authors of event generators often 
choose their own conventions in encoding the details 
of hard process such as qq — > ngZ/"f',Z/"f — > A' + /J 
into the event record. 

The NLO solution for PHOTOS would therefore be 
feasible with some universal, standard event record, 
nonetheless difficult due to practical issues of inter- 
facing. However, the NLO precision in PHOTOS for 
today and tomorrow experiments is most likely not re- 
quired. For the time being the problem remain rather 
academic. That is why we recommend the users to 
remain with the program version as distributed since 
October 2005. Only the users interested in special 
tests and in particular incorporation of decay channel 
dependent electromagnetic formfactors are advocated 
to incorporate further optional changes. In such a case 
even larger care on event record content is required. 

In ref [14], we presented similar modifications in 
the PHOTOS kernel for the decay of B mesons into 
a pair of scalars. The implementation of the ex- 
act (scalar-QED only) kernel brings a minuscule im- 
provement in the agreement between PHOTOS and the 
reference exact simulation of SANC lfl6l . In this case 
both: SANC and PHOTOS are used to simulate single 
photon emission (There exists no reference simula- 
tion with which the multi-photon version of PHOTOS 
could be compared.). 

The way to play with form-factors, originating 
from phenomenological models, and fits to the exper- 
imental data is open now. 



4. PHOTOS debugging tools 

During the years of PHOTOS development, various 
software-related problems needed to be faced by its 
authors. The majority of the workload (and actual 
lines of code) needed to be devoted to the treatment of 
the data stored in a "standard" HEPEVT event record. 

In an ideal situation, if all rules of the HEPEVT stan- 
dard definition were respected, PHOTOS would eas- 
ily identify the branching points (particle decays with 
charged products) in the decay tree, extract the re- 
quired data, then eventually append the generated 
photons (if any) as additional decay products. It 
would also traverse the decay tree to identify all pos- 



sible places where QED corrections might need to 
be generated. Initially, the algorithms employed in 
PHOTOS assumed that the data structure is consistent, 
with all pointers (to "mothers" and "daughters") set 
up correctly. An acyclic tree of 1 — * n (or exception- 
ally 2 — > ri) processes could have easily been navi- 
gated using standard algorithms. 

However, the rigidness of the HEPEVT standard, and 
the lack of possibility of extending it in a consis- 
tent way forced the authors of event generators to 
overload the meaning of the elements of the HEPEVT 
data structure. Certainly, one could consider hav- 
ing bi-directional relations (i.e. mothers pointing to 
daughters, and daughters pointing to mothers" as re- 
dundancy, and the place where additional informa- 
tion (such as spin or colour flow) may be stored in- 
stead. The meaning of the pointers become generator- 
specific and the navigation in such data structure 
could not be performed by a generic algorithm. Pan- 
dora's box of event-record problems has been opened, 
hurting mainly the coordinators of the large experi- 
mental simulation chains. 

The pointers in the HEPEVT structure were not 
the only element that became non-standard. Due 
to evolving needs of physics models (such as big- 
ger number of particles being simulated, or preci- 
sion), HEPEVT data structure has been modified to 
store single- or double-precision data, with various 
array sizes. Dubious matching of the HEPEVT layout 
between simulation blocks became yet more compli- 
cated. 

To alleviate the problem with varying precision and 
layout of HEPEVT, PHOTOS has been equipped with 
a set of debugging and data-interpretation facilities. 
Firstly, it was modified to work on a local copy of 
the event record (the layout of which followed the 
"well-behaved" HEPEVT standard), and have a set of 
functions that would transfer the data between what- 
ever external variant of HEPEVT data structure, and the 
internal storage. Secondly, a set of sanity-checking 
and pointer-reconstructing procedures were applied 
during the transfer of the data between the external 
event record and the internal one. Finally, a debug- 
ging function PHLUPA was provided. It prints out the 
data as interpreted/modified by PHOTOS routines, at 
different steps of event construction. 

Because of the specifics of the PHOTOS algorithm, 
namely massive search and modification of the com- 
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plete event tree, PHOTOS itself has became a debug- 
ging tool for large simulation chains in the experimen- 
tal collaborations. Strengthened with its debugging 
tools, it helped to identify many problems related to 
event grammatic. 

This class of problems will remain with C++ appli- 
cations as well. The origin of the difficulty is physics. 
On theoretical side authors of main Monte Carlo pro- 
grams explore different options on storing informa- 
tion on the process generation, on the other, experi- 
mental users prefer to rely only on the information, 
which can be obtained from the measurement. 

To be able to explore, in such conditions, the 
potential of the PHOTOS algorithm, and to make 
the debugging of the event record data easier, a 
tool: MC-TESTER Q originally developed for tests of 
TAUOLA was adopted. It performs comparison tests 
of distributions of invariant masses produced by two, 
or more, (versions of) event generators. In semi- 
automatic way (thus eliminating the risk of acciden- 
tal programmatic error) it extract the data from event 
records filled by an event generators: it identifies the 
decay modes of a given particle, and for every mode it 
builds the histograms of invariant masses of all com- 
binations of decay products. 

5. Challenges of event record 

On reading the paper, one could get an impres- 
sion that the communication between the modules 
of the simulation tree is a challenging, yet standard, 
goal, which could be realized in any modern, or even 
not so modern software environment. Our discus- 
sion in the previous chapters pointed to possible con- 
straints and requirements in the organization of such 
data structure imposed by relatively modest applica- 
tion PHOTOS, if its precision requirements would need 
to be increased beyond certain level. 

One can ask the question why the C++ implementa- 
tion of PHOTOS did not meet so far as much attention 
as its seemingly obsolete FORTRAN version. On a 
first sight the answer is simple: there was until re- 
cently no commonly agreed standard for the event 
record data structure in C++ accepted by the domi- 
nant part of the community. Recently, it seems that 
the HepMC ifTTl structure is gaining popularity in the 
LHC applications. 

From the past experience of HEPEVT standard, one 



could postulate that a viable event record could be 
seen as a system of parallel trees, the nodes of each 
being bound with correspondence relations spanning 
across the layers. Each individual tree should be easy 
to investigate or modify by program such as PHOTOS, 
or the detector simulation software. 

For MC-TESTER to be effective in comparing the re- 
sults generated by various Monte Carlo simulators, 
it was essential to provide access to various stan- 
dards (or: flavors) of event record data structures. 
Similarly to PHOTOS, MC-TESTER performs exhaustive 
search and data extraction from event record, it how- 
ever doesn't need to modify the contents of the event 
record. Typical data that needed to be extracted was 
the mother-daughter relationships (including finding 
out the non-decaying final-state particles in the cas- 
cade decay), and determining the properties (four- 
momenta, and type) for involved particles. To sepa- 
rate MC-TESTER from problems related to event record 
processing, the HEPEventLib abstraction layer was 
created 

At the technical layer, HEPEventLib does nothing 
more than interpretation of the data stored in various 
event record standards, and providing these data in 
a consistent form to the main program, hence hiding 
all dependencies and data- translation operations. The 
data is provided by abstract object representing a par- 
ticle, a list of particles and an event. As no mod- 
ification is performed on the structure of the event 
record, the properties visible in a "particle" object 
may be mapped directly to the corresponding data in 
the underlying event record. Thanks to that feature, 
particle's properties (such as four-momenta, but not 
the mother-daughter attributions) might even be mod- 
ified from within HEPEventLib's abstract view, and 
the changes would be propagated to the actual event 
record in a consistent way. 

At the bottom of the HEPEventLib there are imple- 
mentations of the HEPEventLib abstraction to con- 
crete event record types, of FORTRAN or C++ appli- 
cations. New "backends" for any (future) event record 
may be implemented as needed - MC-TESTER will au- 
tomatically profit from the new standard with no need 
of adding a single line of code in it. 
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6. Summary 

In the present talk we have reviewed the basic 
properties of theoretical (QED) and software envi- 
ronment which is at foundation of design and per- 
formance of PHOTOS Monte Carlo program for sim- 
ulation of QED bremsstrahlung corrections in decays 
of resonances and particles encoded in different type 
of event records used in simulations of High Energy 
physics. Special emphasis was on the points impor- 
tant for evaluation of systematic errors. 
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